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ABSTRACT 


A one-dimensional boundary layer model is used in conjunction 
with satellite derived infrared surface temperatures to deduce 
values of moisture availability, thermal inertia, heat and evaporative 
fluxes. The Penn State satellite image display system, a sophisticated 
image display facility, is used to remotely sense these various 
parameters for three cases: St. Louis, MO; The Land Between the Lakes, 

KY; and Clarksville, TN. 

The urban centers display the maximum daytime surface 
temperatures which correspond to the minimum values of moisture 
availability. The urban center of St. Louis and the bodies of water 
display the maximum nighttime surface temperatures which correspond 
to the maximum thermal inertia values. It is shov-n that moisture 
availability and thermal inertia are very much responsible for the 
formation of important temperature variations over the urban-rural 
complex. 

Maximum heat flux and minimum evaporative heat flux values are 
found in the urban centers which correspond closely to the sites of 
maximum daytime temperatures. The opposite is found for the rural 
areas. Therefore, a large fraction of the variability in heat and 
evaporative fluxes is caused by variability in land use and in 
rates of evapotranspiration. 

Satellite determined measurements of heat and evaporative 
fluxes and visible channel albedos are compared to corresponding 
values made by aircraft. In the case of surface heat flux the 
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agreemenC was poor; with surface evaporative flux there was some 
agreement; while the albedos agreed very closely. 

The Willis and Deardorff method of estimating dispersion 
parameters is used to infer the effects of surface heat flux on 
the downwind plume spread and concentration of contaminants. The 
ground>level point concentrations are found to be largest for the 
smallest heat fluxes; however, the maximum ground-level point 
concentrations are found to be Independent of heat flux. In view 
of the large variability of heat flux with land use. it is Important 
to consider the land use in the downwind direction of a power plant. 
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1.0 INTRODUCTION 


1.1 StatcMDt of th« Problom 

In th« past f«w years, many studies have been conducted to 

directly measure heat and evaporative fluxes, and substrate 

parameters, such as moisture availability and thermal inertia. 

Spittlehouse and Black (1979) used Bowen ratios and eddy correlation 

measurements to determine the evapotransplration from forests. The 

major problem was that the vertical temperature and humidity 

gradients above the forest canopy were very small, generally less 
-1 -3 -1 

than 0.1*Cm and 0.1 gm m , respectively. These small gradients 
were a result of the high. degree of turbulent mixing generated by 
the large roughness of forests. Small errors in these measurements 
resulted In large errors in the determination of evapotransplration. 
Other studies, including Hicks et al. (1975), also found it difficult 
to measure forest evapotransplration. 

On 26 April 1978, the National Aeronautics and Space Administration 
(NASA) launched a small Applications Explorer satellite as part of 
the Heat Capacity Happing Mission (HCMM). The satellite carried a 
tm-channel radiometer (0.5*1.1 ym and lO.S^lI.S ym) in a sun- 
synchronous orbit. The satellite had a resolution of .6 km and 
collected radiometric ground temperature data at approximately 0230 
and 1330 LST which was close to the times of maximum and minimum 
surface temperatures. The satellite capabilities are discussed in 
detail by Price (1979). Dodd (1979), Kocin (1979), and Carlson et al. 
(1980) used data from the HCMM mission in conjunction with a one- 


dloMMlonAl boundary layar nodal davalopad by Carlton and Boland (1978) 
to infar valuaa for noiatura availability* tharmal Inartia* haat 
and avaporatlva fluxaa. Tha aathod uaad by thaaa workara to 
datamlna such paranatara fron ramotaly sansad data raqulrad a 
consldarabla amount of tlna for procassing tha convantlonal batch 
computar data. 

Once the problem of obtaining rali.*bla surface haat flux data 
over a particular area has bean solved, one application Is to study 
the effect of a variable surface heat flux on plume spread and 
concentration. This can have a profound effect on the downwind 
spread and concentration of the plume. 

1.2 Purpose of the Thesis 

This research applies high-resolution surface temperatures 
obtained by the use of a sophisticated satellite image display system 
in conjunction with a one-dimensional boundary layer model to infer 
values of moisture availability, thermal Inertia, heat and evaporative 
fluxes. The Penn State satellite image data manipulation system, 
discussed in detail in section 2.0, involves a faster, more reliable 
method for determining the surface and substrate parameters than 
that used in earlier HCMM efforts. 

Three case studies, St. Louis, tff); The Land Between the Lakes, 

KY; and Clarksville, TN, are presented in sections 3.0 and 4.0. The 
geographic locations of these sites are displayed in Figure 1. The 
St. Louis case study was conducted for the purpose of investigating 
the surface energy balance and substrate parameters over a large urban 



Figure 1. The geographic locations of Che three case studies: 
SC. Louis, MO; The Land Between the Lakes, KY; and 
Clarksville, TN. 
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ar««. The results ere coopered with pest urben heet Islend reseerch. 
The Lend Between the Lekes end Clerksvllle cese studies were 
conducted es e pert of the Envlronmentel Protection Agency's STATE, 
"Sulfur Trensport end Trensforoetlon In the Environment", project. 
Albedo, heet end eveporetlve flux i&eesureaents obtslned by slrcreft 
ere compered with the setelllte meesured velues. 

The surfece heet flux velues for the three cese studies ere 
exemined to see how the surfece heet flux verles over different 
terreln. The Willis end Deerdorff method of estlmetlng dispersion 
peremeters Is then used to study the effect of this verleble surfece 
heet flux on plume spreed and concentretlon. Most Importently, the 
varletlon of the ground-lrvel point concentrations with differing 
heat fluxes Is examined. This Is presented In section 5.0. 
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2.0 THE PENN STATE SATELLITE IMAGE DATA MANIPULATION SYSTEM 

2.1 Past Method of Analysis 

The sequence of computer programs used by Augustine (1978) , 

Dodd (1979), and Kocln (1979) to extract the raw data from the 
computer tapes supplied by NASA and to grid the location of the 
working area has been superseded by a sophisticated image display 
system. This system, described in detail in section 2.2, saves a 
tremendous amount of time and increases the accuracy of the final 
product. 

Augustine (1978) first used the program RECSTRA4 to extract a 
designated portion of a HCMM satellite picture from the raw data 
tape to a 512 by 512 area storage tape. The proper scan line and 
pixel from which the subset was to begin had to be specified by 
the user. This process was then repeated using the TGRFLOAS 
program to subset a 130 by 130 disk file from the 512 by 512 tape 
file. A histogram was also simultaneously produced. The histogram 
was used to assign characters equivalent to the digital data on disk. 
A character map was then produced which gave the user his first view 
of the area subsetted. If this area was not exactly what the user 
wanted, he would have to go through this time-consuming process 
again. Then the user would have to adjust the character assignments 
to obtain the best possible character map of the area desired. 
Identifiable physical features would have to be found from the 
character map to be used in conjunction with known reference points. 
This would involve the tedious process of counting the number of 
pixels In the X and Y direction from the character map. 
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A calibration equation for the HCMM data was used to convert 
the digital count or DN values (0-255) into temperatures (260-340°K) . 
A water vapor corre:tion suggested by Cogan and Uilland (1974) was 
also applied to the calibration equation; however, the water vapor 
correction is now calculated through the program called NASRAD 
which is discussed in section 2.3. 

The final analysis was then displayed through the use of a 
CALCOMP plotter to produce, for example, an Isotherm analysis. The 
false-color-enhanced Image created by a RAMTEK color display monitor 
was also used to display the final analysis. 

This whole process has now been updated, so that the task of 
extracting, mapping, and contouring a desired area takes only one 
or two days as opposed to two or three weeks. The method of deriving 
the maps of SKJisture availability, thermal Inertia, heat and 
evaporative fluxes is still used and this method is discussed in 
detail by Dodd (1979). 

2.2 Data Display and Analysis 

The purpose of the Penn State Satellite Image Data Display 
System is to enable the user to extract a workable area from a 
TIROS-N or HCMM satellite image and to transmit this area to an 
IBM 370 BAT (Batch and Terminal) file. A three step process is 
used to accomplish this. First, the data are stripped from the 
master data tape obtained from NASA and written to a scratch tape. 
HCMM data are written to the scratch tape as it appears on the 
master tape. TIROS-N data are calibrated and converted from thermal 
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energy units to absolute temperature, and are normalized on a scale 
of 0 to 255 where 0 is equal to -30*C and 255 is equal to 70*C. 

The scratch tape image is then displayed on a video monitor, and a 
smaller, workable area is selected for transmission. Finally, the 
data are transmitted to an active file on the IBM 370 via Penn 
State's RJE (Remote Job Entry) system. 

The following equipment is used to extract and transmit an 
image to an IBM 370 BAT file: 

1) PDP 11/10 minicomputer operating under RT-11 
version 3.2 

2) PDP 11/34 minicomputer operating under RSX-*11M 
version 3.2 

3) RK05 and RK07 disk drives 

4) Digi data nine-track tape drive 

5) Grlnnell system (processor unit and T.V. monitor) 

6) DR 11-B interface (direct measuring access) for 
the Grinnell processor 

7) Grlnnell monitor control box containing two 10-turn pots, 
one button, and one toggle switch 

8) Modem (Dataphone 300/1200). 

Once the master tape has been stripped, the communications file 
program called COMFIL is used to subset and transmit a workable area. 

A workable area usually encompasses approximately 4225 square 
kilometers. A series of commands are used to accomplish these casks, 
and they are described below: 

1) HELP — Used to obtain a description of and Instructions 
for running any of the following commands. 
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2) TAPE DISPLAY — Displays tha scratch Capa imaga 
(or portion of It) on tha Grinnall monitor. 

3) GRAY SCALE — Allows tha usar to adjust tha 
resolution of the screen Image. 

4) COORDINATES — Used to select a portion of the 
screen Image to be enlarged or stored on disk. 

5) READBACK — Types the tape pixel coordinates and 
Intensity of selected points on the video Image. 

6) SAVE TAPE — Saves that portion of the tape Image 
specified by the COORDINATES consnand. 

7) HISTOGRAM — - Constructs a histogram of pixel 
Intensities for the Image file created by the 
SAVE TAPE or CLOUD command. 

8) CLOUD — Removes clouds from the Image saved by 
the SAVE TAPE command. 

9) DISK DISPLAY — Displays on the T.V. monitor the 
disk file Image created by the SAVE TAPE or CLOUD 
command. 

10) CONVERT RJE — Converts a disk file Image Into a 
form suitable for transmission to Penn State's RJE. 

11) RECOVER — Used to recover from RJE transmission 
errors. 


A schematic representaclon of this analysis procedure Is 
depicted in Figure 2. The Image first displayed by the TAPE DISPLAY 
command is 2048 by 2048 pixels and an example of such an image Is 
shown in Figure 3. This Is a day Infrared HCMM image for 22 August 
1978 taken at 1358 LST. The area extends from Indiana and Ohio to 
Tennessee. The black areas indicate low temperatures; whereas, 

Che white areas indicate high temperatures. The blackest regions in 
Figure 3 reveal cirrus clouds to the north and cumulus clouds to 
the south. The urban heat islands of Indianapolis, Indiana; 



Figure 2. Schemacic representation of the Penn State 
satellite Image data system. 







Figure 3. The 2048 by 2048 pixel laege for 22 August 

1978 et 1338 LST. The derkeet ereee Indlcece 
low tei^ereturee end the lightest ereee 
indicete high teaperetures. 


Figure 4. The 512 by 512 pixel Imegc choeen in Figure 3 
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CiaeloMCl, Dayton« and Colui^u«i Ohio art racotniiabla at tha 
whitatt araaa in tha norcharn portion of cha iaaga. Tha Land Batwaan 
tha Lakaa, Kantucky, and Clarkavilla, Tannaaaoa* tha caaa atudiaa 
daaeribad in aaetion 4.0, ara viaibla in tha boxad araa of tha photo. 

Tha COORDINATES prograa ia usad to axtraet tha aubaattad iaaga. 

Tha haight and width in pixala ara apacifiad and an appropriataly 
acalad box containing croaahaira ia aupariapoaad onto tha I'ccaan iaaga. 
Tha croaahaira aupariapoaad onto Flgura 3 ara 512 by 312 pixala. Tha 
box ia aovad by turning tha X and Y knoba on tha control box. TAPE 
DISPLAY ia than uaad to diaplay thia aaallar araa. Figura 4 diaplaya 
tha araa choaan to ba anlargad in Flgura 3. Thia aaaM proc^^aa ia uaad 
again to chooaa tha final, workabla araa. Tha box aupariaq>oaad onto 
Figura 4 ia 128 c;* 128 pixala and Figura 5 ahowa this final araa. At 
any tioia, tha GRAY SCALE coonand can ba uaad to adjuat tha raaolution 
of tha acraan inaga. 

Thia inaga ia now aavad on diak by uaa of tha SAVE TAPE conaMnd. 
HISTOGRAM ia chan uaad to conatruct a hiatogran of pixal intanaitiaa 
which providaa infomation about tha tanparatura diacributiop in 
cha working araa. Any abnomally cold apikaa on 'ha hiatogran 
indicata chara ara clouda in tha final 128 by X2i^ avaa; tharafora, 
it is nacassary to uaa Cha CLOUD REMOVAL coonand. A chin cirrus 
cloud is visible in cha upper rlght~hand corner of Figura 5. Tha 
low and high indices of cha cloud, obtained fron Cha hlstogra.^^, 
along with a default value ara antarad with Cha CLOUD REMOVAL coonand. 
If the actual radius of a cloud is lass chan or aqua} Co the sracifiad 
value, Che cloud will be pacched over by tha intensity values of tha 


Figure 5. The 128 by 128 pixel image chosen in Figure 4 


Figure 6. The 128 by 128 pixel image chosen in Figure 4 
with the cirrus cloud in the northeast portion 
of the image removed. 
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plxsls on its psriphery. If ths cloud radius is larger then ths 
supplied value, ths cloud will bs' patched over by ths default value, 
which is usually vary black (intensity value 0) or very whits 
(intensity value 255). This process is repeated until the cloud, 
such as the one in Figure 5, is renoved. This cloudless image, 
which is shown in Figure 6, is then resaved on disk using the SAVE 
TAPE coonand. 

The latitude and longitude (navigational coordinates) of the 
corner points are found via the READBACK command. READBACK types 
the tape pixel co:>rdlnates and intensity of selected points on the 
T.V. image. A set of crosshairs are displayed on the screen and 
they are positioned over prominent landmarks by turning the X and Y 
knobs on the control box. The tape coordinates with respect to 
the upper left-hand corner of the screen and the intensity of the pixel 
beneath the crosshairs are found by pushing the button on the control 
box. The READBACK reference points are matched with points of known 
latitude and longitude (reference points) . The utility program 
called LONLAT, described by Augustine (1978), derives linear 
regression equations which describe the variation of latitude and 
longitude along the four borders of the 130 by 130 matrix. LONLAT 
then applies these equations to find the navigational coordinates 
of the corner points. The many Inlets and peninsulas of the Land 
Between the Lakes are well represented, so the reference points are 
easily located and LONLAT la easily applied. 

The 128 by 128 matrix is then converted into a form suitable 
for transmission by use of the CONVERT RJE command, and finally 
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transmitted to an actlva flla on tha IBM 370 via Pann Stata'a BJB 
ayatam. This flla la tranamlttad as a 130 by 130 matrix to ba 
co^>atlbla with LONLAT and othar utility programs. 

This vhola procass la than rapaatad for a night Infrarad Imaga. 

Tha utility program callad REGGIE» dlscussad In datall by Dodd (1979) > 
la usad to extract a region common to both tha day and night working 
areas. REGGIE utilizes a calibration equation for the HC31M data 
which converts the DN values (0>255) Into temperatures (260-340*K). 

Also Included in this calibration equation Is a water vapor correction 
that is calculated through the program called MASRAD. NASRAD coiq>utes 
the temperature difference between the observed ground surface 
temperature and the temperature measured by the radiometer In the 
11 v> atmospheric window. The attenuation of energy Is computed by 
applying the radiative transfer equation to an atmosphere consisting 
of layers for which the average temperature, pressure, and humidity 
are specified. 

Before the final analysis can be made, the output from REGGIE 
must be coupled with output from the model which produces temperature 
as a function of moisture availability and thermal Inertia. The 
model is discussed briefly in section 2.3. 

2.3 The Model 

Boland (1977), Dodd (1979), Kocln (1979), and Carlson et al. (1980) 
have thoroughly documented Carlson and Boland's technique of inferring 
moisture availability, thermal inertia, heat and evaporative fluxes. 

This method is briefly discussed below. 
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Th« oiM*dlffl«n«lonal boundary layar nodal waa davalopad to 
Invaatlgate the liqtortance of varloua terrain and atnospherlc 
varlablea on the aurface tenperature and heat flux. Honin<>Obukhov 
theory la used vlth an implicit K-type parameterization for the eddy 
fluxes in the surface layer. Net Radiation (R^) is the primary 
forcing with no advection or anthropogenic heat sources included: 

R • H + E + G . (1) 

n o o 0 

The surface heat flux (H^) is expressed as a function of an 
eddy and molecular diffuslvity, and a vertical potential temperature 
gradient. The water vapor flux (E^) Is expressed as a function of 
an eddy and molecular dlffuslvlty for water vapor, and a vertical 
mixing ratio gradient. The ground flux (G^) is parameterized in 
terms of the soil thermal conductivity and the temperature gradient 
Immediately below the soll-alr Interface. 

The net radiation Is calculated from the radiative fluxes at 
the surface by: 

*n • ® 

where S Is the solar flux absorbed at the ground, Is the downward 
terrestrial flux, and Is the upward terrestrial flux. Augustine 
(1978) gives a detailed explanation of the radiative transfer model. 

The basic structure of the model which Includes a substrate 
layer, transition layer, surface layer, and mixing layer Is depicted 
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in Figure 7. Input parameters for the laodel include the temperature 
of the lowest level in the ground ('v 1^5 m), roughness length, albedo, 
surface wind speed and direction, preclpitable water, ground 
amlasivity, and turbidity of the atmosphere. A sounding of the 
atmospheric temperature and dewpoint extending from the surface to 
approximately 600 mb la an Important input parameter. An average 
sounding of the day desired and the day before is used as initial 
conditions. 

The mixed layer height (h) is calculated through a formulation 
developed by Tennekes (1973) in %fhich h grows throughout the day as 
a consequence of the surface heat flux (H^) from below and entrainment 
of air with a higher potential temperature from above. The heat flux 
becomes dependent on the temperature lapse rate near the ground at 
night when the upward heat flux disappears. The model then uses a 
method developed by Blackadar (1979) where the nighttime turbulence 
is calculated as a function of the Bulk Richardson Number in the 
surface layer. This formulation is used to simulate the boundary 
layer at night during which turbulent episodes occur. 

The model is initialized at 0600 LSI and allowed to proceed for 
a day. The model is reinitialized, after each cycle, with incremented 
values of moisture availability (.05 - 1.0) and thermal inertia 
(.005 - 0.1) until sixteen cycles are completed. 

Moisture availability (M) is the fraction of maximum possible 
evapotransplration for a saturated surface. Tanner and Felton (1960) 
define M as 

E 

o,pot 


( 3 ) 
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where Is the surface evaporation rate and Is the potential 

evaporation rate or the rate of evaporation from a saturated surface 
at temperature Ik)lsture availability takes into account the 
reduction in the evaporation efficiency due to subsaturation of the 
surface air. The range of M Is near zero for completely dry surfaces 
to 1.0 for saturated terrain. M Is initialized In the model and 
held constant throughout the day. This is in agreement with Nappo 
(1975) who found that the specification of a diurnal range for M 
offers little advantage In modeling the surface fluxes. 

Idso et al. (1976) maintain that Internal and external factors 
must be taken Into account when one looks at the diurnal surface 
temperature variation of any substrate material. Thermal conductivity 
(X), density (p), and specific heat (c) of the substrate are internal 
factors which are used to define thermal Inertia (F) : 

P • (Xpc)^'^^ (4) 

-2 -1 - 1 ^'’ -2 -1 1/2 

which has the units of cal cm K s ' " or Wm K s ^ . The 

latter units are called thermal inertia units (TIU) where approximately 

“2 -1 - 1/2 

400 TIU are equivalent to .01 cal cm K s . Solar radiation, 
air temperature, atmospheric precipitable water content, wind, etc. 
are some of the external parameters which affect the amplitude of the 
diurnal surface temperature wave. 

As will be shown, the thermal inertia is most responsible for 
shaping the nighttime temperature pattern; whereas, the moisture 
availability has a greater effect on the temperatures during the day. 


21 


3.0 THE ST. LOUIS URBAN HEAT ISLAND 

3.1 Past Urban Heat Island Research 

Man has altered his environment in many ways with urbanization 
and industrialization being the best examples. The fact that urban 
areas are warmer than their rural surroundings has been known for a 
very long time. Luke Howard first published evidence of the urban 
heat island in 1833 when he showed that air temperatures measured 
in London were about 1.1°C higher than those in outlying areas. In 
1927, Wilhelm Schmidt was the first to use an automobile to obtain 
a temperature analysis of a city. Since then, many studies have 
been conducted to determine the climatic impact of altering the 
terrain. 

There are many reasons why cities possess a different climate 
than their environs. Landsberg (1979) maintains that the radical 
change in surface characteristics is most important in creating the 
heat island. Cities are made up mostly of brick, concrete, and 
asphalt which are believed to have a high heat conductivity and 
capacity. The heat stored during the day is slowly released at 
night causing warmer surface temperatures in the city. The anthro- 
pogenic heat produced from combustion (space heating, manufacture, 
transportation, lighting) can contribute 13 to 30% of the total heat 
budget. In fact, in Manhattan this heat source can contribute two 
and a half times the solar input in the winter. The tall buildings 
of the city increase the chance of multiple reflections of radiation. 
There is also a reduction in albedo, evaporative heat flux, and 


ventilation. 
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Duckworth and Sandberg (1954) atudled three California cltiea 
of various sizes. The temperature was found to increase from rural 
to urban centers in direct proportion to the structure density of 
the city. Chandler (1965) conducted a 30 year study of the climate 
of London. He found that the city had more days with exceptionally 
high temperatures and fewer with low temperatures cosq>ar«d to the 
rural surroundings. 

Mitchell (1961) gives three reasons why urban heat islands are 
not topocllmatic anomalies that would exist even in the absence of a 
city. First, cities of completely different topography possess heat 
islands. Second, the heat Island Is comparatively weaker on Sundays 
when Industry and other urban activities are at a minimum. Third, 
Mitchell found that the intensity of the urban heat island increases 
with Increased size and population of a city. 

In 1971, an extensive five-year investigation called METROMEX 
was undertaken to see how St. Louis affects summer weather. White 
et al. (1978) found an extreme difference of 1S*X in maximum infrared 
surface temperatures between industrial and undeveloped land. Solar 
noontime albedo values were found to vary from 15-172 for rural land 
In contrast to 12-132 for urban land. Measurements showed that the 
energy loss by emitted radiation is Inversely related to the amount 
of vegetation in the St. Louis area during mid-August. Changnon 
(1978) found that within 40 km of the city in the summer, increases 
in the frequency of occurrence of various thunderstorm characteristics 
("V 4-10 to 1152), hailstorms (+3 to ■♦•3302), various heavy rainfall 
characteristics (+35 to +1002), and strong gusts (+90 to +1002) occurred. 
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Braham and Wllaon (1978) found from radar acho Maaurananta that 
St. Loula axarta a aubatantlal Influanca on tha halghta of convactlva 
clouda. Thara la roughly a twofold probability of finding tha tallaat 
acho ovar and downwind of tha city In tha tvmmr during tha nidday 
and lata avanlng houra. 

It haa only baen during tha laat alx yaara that aatalllta 
radiation data hava baan uaad to atudy tha urban haat la land. Carlaon 
at al. (1977) uaad data from tha MOAA-3 VHRR (Vary High Raaolutlon 
Radlomatar) to atudy tha Infrarad radlatlva (blackbody) aurfaca 
tanparatura dlatrlbution for tha Loa Angalaa araa. Tha hlghaat 
morning tamparaturaa wara found ovar tha Induatrlal zona, and In tha 
avanlng tha hlghaat tamparaturaa occurrad ovar tha cantral bualnaaa 
diatrict and raaidantlal araaa. Mataon at al. (1978) uaad tha HOAA-5 
aatalllta to datact mora than fifty urban haat lalanda in tha mid- 
waatarn and northaaatam Unltad Stataa on a cloud-fraa aummar avanlng. 
Tha maximum urban-rural tamparatura diffarancaa rangad from 2.6*C to 
6.5*C. Carlaon at al. (1980) uaad HCMM Infrarad aatalllta tamparatura 
maaauramanta In conjunction with a ona-dlmanalonal boundary layar 
model to infer tha dlatrlbution of aurfaca haat and evaporative 
fluxaa, ground molatura availability, and thanaal Inertia. A 
reduction of evaporation and molatura availability along with an 
increaaa of aanaibla haat flux waa found ovar urbanized araaa and 
ovar cropped araaa with low vegetative cover. The thermal inertia 
distribution waa ill-defined and its variation between urban and 
rural areas was quite small. 
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This h«t b««n • v«ry brief review of peer urban Kent ielend 
reeeerch. A detailed review of recent urban ellaatology atudlea 
la given by Oke (1979a). 

3.2 Background Auguat 23-24, 1978 

St. Louie, the largest city of Missouri, is situated on the 
west bank of the Mississippi River approxiaately 16 km south of 
its Junction with the Missouri River. St. Louis is one of the chief 
cities and major industrial centers of the midwest. The city has a 
population of approximately 623,000 making it the eighteenth largest 
city in the United Stetes. The city is spread along 31 km of the 
Mississippi River and rises gently to the tfest to an elevation of 
142 meters. St. Louis county, which extends all along the city's 
western perimeter, has more rolling terrain. Except for industrial 
parks, most of St. Louis county is residential, with the western 
portions being somewhat rural. The other counties are more rural 
and forested. In contrast, most of the St. Louis sMtropolitai. area 
on the Illinois side of the Mississippi River is industrial. The 
urban areas there are partly flood plain; however, there are bluffs 
in some areas. The Isrgest urban center on the eastern bank of the 
Mississippi River is East St. Louis vdiich has a population of 
approximately 70,000. 

The weather in St. Louis on 23-24 August 1978 was dominated by 
a large ridge of high pressure centered over West Virginia. Clear 
skies and tegq>eratures around 35 *C and nighttime temperatures around 
18*C were recorded at the airport. Precipitation had not occurred 
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ov«r tiM 8e. Lottit «rM for at laaat a waek; tharafora* tha gtcmd 
can ba axpaetad to hava baan qulta dry. 

Tba aatallita tOMgat wara takan twalva houro apart on 23 August 
1978 at approxiaataly 0230 L8T and 1330 L8T. Tba nodal output uaa 
for 1330 LST on 23 August and 0230 L8T on 24 August; thus. It was 
assunad that tha nlghttlaw satalllta inagt oecurrad 24 hours aftar 
Its satalllta orbit tltM. It Is ballavad that this ravarsal of 
orbit tlnas did not affact substantially tha final rasults. Carlton 
at al. (1980) assunad ttatlonarlty la tha tanparatura rasponsa 
batwaan two days as a basis f )r ravarsing tha ordar of tha day/nlght 
taquanca In tha nodal. 

Sons of tha Initial paranatars usad for tha St. Louis nodal run 
ara llscad In Tabla I. Tha 1200 QfT soundings for 23 and 24 August 
1978 ara shown In Flgura 8. Tha avaraga of thasc two soundings, 
also saan In Flgura 8, was usad as input to tha boundary layar nodal. 

Tha base map along with key areas chat hava baan coded for easy 
refaranca Is displayed in Flgura 9. 

3.2.1 August 23 Daytlsa Tanparacures 

Tba daytlna surface blackbody taaq>araturas displayed In Flgura 10 
show a wide spread of values ranging from 45*C to 29*C. As expected, 
tha urban areas show tha highest caaparaturas. Surface temparaturas 
greater chan 43*C ara recorded in downtown Sc. Louis (D) with tvo 
maxima of 45*C centered along tha interstate highway. Granite City (C), 
a highly industrialized area, exhibits a maximum of 45*C and an area 
of 45*C is also cantered over the urbanized area just south of Jennings. 
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T«bl« I. Initial Modal Input Paranatara fov St. Loula 


Ground Taa^ratura at tha 290*K 

Basa of tha Layar 

Roughnaat Langth 25 ca 

Praclpltabla Watar 2.54 ca 

Albado .19 

Ealatlvlty 1.0 

-1 

Wind Spaad at a ^ 3 a 3.6 asac 

Dlractlon 210* 



Figure 8. The 1200 GMT Sc. Louis soundings for 23 and 24 August 
1978, along with the average of the two soundings. 


Figure 9. The base map used for the St. Louis case study along 
-rith the following key areas coded for easy reference 

D Downtown St. Louis 

E East St. Louis 

P Forest Park 

H Horseshoe Lake 

G Granite City 

C Crop and Pasture Land 
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Th« highly industrialized end urbanized area of East St. Louis (E) 
shows a temperature maximum of A3*C. An area of 41*C is also seen 
just south of Horseshoe Lake over Fairmont City, an industrial and 
residential area. 

The physical and thermal properties of the city, along with the 
heat derived from induatrlal space heatings contribute to the high 
temperatures. Areas of maximum temperature outside the urban complex 
may be associated with industrial parks. Nkemdlrim (1976) maintains 
that these high temperatures are partly due to the enhanced transfer 
of heat into the atmosphere through chimneys, the lower aggregate 
specific heat of the industrial park fabric, heavy truck traffic, 
and the enhanced horizontal mixing of air currents due to low 
density developments. 

A 6"C variation of temperature is exhibited south of Horseshoe 
Lake (H) where the land use changes from an urban to a rural fabric. 
The lower temperatures are recorded over cropland, pasture, and 
wetland areas. Moreover, the crop and pasture lands (C) north of 
the city show temperatures ranging from 10*C to 12*C lower than 
the downtown area. The lowest temperatures are recorded along the 
Mississippi River «diere three minima of 29”C are located. A minimum 
of 31 is also located over Horseshoe Lake. 


3.2.2 August 24 Nighttime Temperatures 


The range of nighttime temperatures displayed in Figure 11 is 


much smaller than the daytime range. Two small areas of 19"C are 
found along the Mississippi River. Areas of 18°C are found over the 
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downtown area, Horsashoa Laka, Granlca City, and Eaat St. Loula. 

The fimall 18*C maximum over the northwest corner of the map is 
associated with a narrow area of conaercial activities, as is the 
small 18*C maximum south of the downtown area. The heat stored 
during the day in these areas is slowly released at night causing 
warmer surface temperatures. Duckworth and Sandberg (1954) point 
out that the heat capacity of urban building materials as compared 
with that of soil to release stored heat energy during the night is 
sufficient to account for well over half of most observed urban 
temperature differences. The rural areas north and east of the 
city display minima of lA^C to 15°C. However, the model does not 
Include anthropogenic heat effects. 

3.2.3 August 23-24 Moisture Availability 

The moisture availability (M) values displayed in Figure 12 
range from .1 over the urban areas to .7 over some rural areas. 
Values of M that approach 1.0 are found over the bodies of water. 

Minimum values of .1 are found over the downtown area. 

Granite City, East St. Louis, and the area south of Jennings. 

This is obviously due to the large areas of concrete and lack of 
vegetation. Small values of M are also found over some cropland 
areas where the vegetation is sparse; whereas, larger values of M 
correspond to more vegetated areas. Forest Park (F), just east 
of the downtown area, displays a value of .3. Larger values of M 
are located here but smoothing of the data reduced the displayed 
value. This is one of the largest municipal parks in the nation. 
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so Its size of 1400 acres should show up with a somewhat higher 
value of M than the areas of concrete surrounding it. The maximum 
value of .8 over Monsenthein Island, northwest of Granite City, is 
due to the fact that a dense canopy of deciduous trees is prevalent 
here. 

Moisture availability is inversely related to the daytime 
temperature. For example, .1 corresponds with the 43"C daytime 
isotherm: similarly, .8 corresponds with the two 29*C isotherms 
over the Mississippi River and the 31“C isotherm over Horseshoe Lake. 

3.2.4 August 23-24 Thermal Inertia 

The values of thermal inertia (P) depicted in Figure 13 range 
-2 -1 - 1/2 

from .015 to .07 cal cm K s . Large values of P are 
associated with the rivers and Horseshoe Lake; however, the values 
of P are not reliable over water. 

The thermal inertia is closely correlated with the nighttime 
tamperatures. Low values of P (.02-. 03) are found north and east of 
the city over rural areas where the ground does not effectively 
store heat during the day and does not conduct it to the surface 
at night. Thus, low nighttime temperatures (14-16*0 correspond 
with these low values of P. Conversely, high values of P (.04-. 05) 
are found over urban areas near downtown and East St. Louis where 


there is a high heat capacity of the surface. It is here that 
the highest nighttime temperatures are found. 


ST. LOUIS, MO I «LJUn 1 


Figure 13. Thermal inertia analysis for 23-24 August 1978 
(cal cm-2 k“1 s-1/2). 
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3.2.5 August 23 Surface Hcet Flux 

The surface heat fluxes (H^) shown in Figure 14 correspond to 

the early afternoon of 23 August. High values of heat flux correspond 

-2 

to high daytime temperatures. Heat flux maxima of 200 Wm are 

found near downtown St. Louis » south of Jennings, and Granite City. 

These isopleths correspond with the 45”C daytime isotherms. East 

_2 

St. Louis also displays a large heat flux of 180 Urn which corres- 
ponds to the 43*C daytime Isotherm. Conversely, low values of heat 

. .ux correspond with low daytime temperatures. Heat flux minima of 
-2 

60 Wm are found over Horseshoe Lake and the Mississippi River. 

These isopleths conform with minimum daytime temperatures of 29*C 
to 31“C. 

This variability of heat flux with land use corroborates the 

findings made by Ching et al. (1978). In their study, sensible 

heat fluxes were calculated by the eddy correlation method. During 

the summer, the heat flux was found to be largest at a site composed 

mainly of paved areas and other dry surfaces. Figure 14 shows that 

the largest heat fluxes in downtown St. Louis and south of Jennings 

are found over interstate highways which traverse northwest frcna 

the city. Ching et al. (1978) also found smaller heat fluxes over 

more vegetated urban sites with active evapotransplration. 

«2 

High values (200 Wm ) also correspond to low values of 

-2 

M (.1): likewise, low H^ values (60 Wm ) correspond to high values 
of M (.8). 


ST. LOUIS, MO 


Figure 14. Surface heac flux analysis for the early afternoon of 
23 August 1978 (Wra“^). 
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3.2.6 August 23 Surface Evaporative Heat Flux 

The values for the evaporative heat flux (E^) shown in Figure IS 
_2 -2 

range from 200 Wm in the urban areas to 445 Wm over the bodies 

of water and the outlying rural areas. The minimum values of E^ 

correspond to the minimum values of M (.1) end the maximum values 
_2 

of H (200 Wm ). Likewise, the maximum E. values correspond to 
o o 

-2 

the maximum M values (.8) and the mini’wm H^ values (60 Wm ), 

o 

-2 

Local maxima such as the 305 Wm isopleth over Forest Park 

reflect the transition from concrete to vegetation. Oke (1979b) 

maintains that the evaporative cooling associated with such surfaces 

as Forest Park may be significant in providing local amelioration 

of heat Island influences. Outlying areas which contain large 

-2 

amounts of vegetation exhibit values that approach 480 Wm . 

3.3 Summary of Results 

The maximum daytime surface temperatures corresponded to the 
minimum values of moisture availability and this was found in the 
center of St. Louis. The maximum nighttime surface temperatures 
corresponded to the maximum thermal Inertia values and this was also 
found In the urban centers. Thus, it was seen that M and P were 
very much responsible for the formation of Important temperature 
variations over the urban-rural complex. 

Maxima of heat flux and minima of evaporation were also found 
In the urban center where the maximum daytime surface temperatures 
were located. The suppression of evaporation on dry surfaces 
allowed a larger portion of the net radiation to be partitioned Into 
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ST. LOUIS, MO 


Figure 15. Surfscs tvapotstivt heat flux for the early afternoon 
of 23 August 1978 (Vhi"^) . 
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••nsiblc hMt* thus causing grsstsr surfscs hasting. Tha opposlta 

was trua for tha rural araas. Mora anargy was partltlonad Into 

latent hast flux for aolst surfacas causing supprassad hasting. 

H and E wars also Inaxtrlcably tlad to tha vagatatlon growth 
o o 

over a particular area. Ralatlvely large hast fluxas and ralatlvaly 
low avaporatlva hast fluxes wars found ovar araas with llttla or 
no vagatatlon. Low hast fluxas and large avaporatlva hast fluxas 
ware found ovar vagatatad araas. Hanca, a larga fraction of tha 
land use variability of hast flux was caused by variability In rates 
of evapotransplratlon, which was caused by the variability of 
transpiring surfaces. 


I 
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4.0 PROJECT STATE RESULTS 
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4.1 Project STATE — The Tennessee Plume Study 

The Environmental Sciences Research Laboratory (ESRL), a 
division of the EPA's Office of Research and Development (ORD), 

Is conducting various research projects designed to provide data 
for the development of regional atmospheric deposition models. 

STATE, an acronym for "Sulfur Transport and Transformation In the 
Environment", was carried out In 1978 primarily for the purpose 
of quantifying the impact of various regional sources on air quality. 

During August 1978, the STATE program focused on the Tennessee 
Valley Authority's (TVA) Cumberland power plant in north-central 
Tennessee. Schiermeier et al. (1979) give the following reasons 
why the Cumberland steam plant was chosen. First, the power plant 
already possessed extensive pollutant monitoring facilities and 
meteorological forecasting data. Also, the site is located within 
the Sulfate Regional Experiment (SURE) monitoring network. The 
plant has a generating capacity of 2600 MW with two 1300 MW boilers 
each connected to its own 305 m stack by three ducts. The three 
ducts are equipped with two electrostatic precipitators per duct 
with sampling ports before and after the precipitators. The coal 
used by the plant is fairly homogeneous and it has an approximately 
3% sulfur content. Finally, Cumberland was chosen because it is a 
base load station which reduces the pollutant emission fluctuations. 
TVA, other government agencies, universities (including Penn State), 
and private organizations conducted the field experiments which were 
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designed Co provide Inforsiaclon on end Insight into the dynsmics of 
plume transport and transformation over long distances. 

4.2 Background August 22-23, 1978 

The Tennessee Plume Study was conducted in northwestern Tennessee 
and southwestern Kentucky. Satellite derived measurements of 
temperature, heat and evaporative fluxes, thermal inertia, and 
moisture availability were made for two adjacent areas: The Land 

Between the Lakes, Kentucky, and Clarksville, Tennessee. 

In Che 1960's, the TVA launched an aad>itious project for a 
national outdoor recreation center located in the Land Between the 
Lakes. The base map used for the Land Between the Lakes study, 
along with key areas that have been coded for easy reference, is 
shown in Figure 16. Kentucky Lake and Lake Barkley are large, 
artificial lakes formed by Che dams of Che Tennessee Valley 
AuChoricy and the U.S. Corps of Engineers. The Land Between the 
Lakes is a heavily forested area with predominantly deciduous trees. 

The base map used for the Clarksville study, along with coded 
areas, is shown in Figure 17. Clarksville, Tennessee, situated 
64 km northwest of Nashville, has a population of approximately 
32,000. Situated along the Cumberland River, the city lies in a 
dark-tobacco, livestock, and farming region. Clarksville is a 
leading marketing center for tobacco products, shoes, and clothes. 

Fort Campbell, 14 km northwest of Clarksville, is on the Tennessee- 
Kentucky border. This U.S. army post, whi.;h has a population of 
approximately 23,000, is the home of the lOlst Airborne Division. 


Figure 16. The base map used for Che Land Between Che Lakea 

case study along with Che following key areas coded 
for easy reference: 

LBL Land Between the Lakes 

KL Kentucky Lake 

JC Jonathan Creek 

BR Blood River 

LB Lake Barkley 

LR Little River 

DC Donaldson Creek 

NSP Numerous Small Ponds 

FCR Fort Campbell Restricted Area 

F Forest and Brushwood 

C Cropland and Pasture 
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Figure 17. The base nap used for the Clarksville case study 
along with the following key areas coded for easy 
reference: 

CR Cuid>erland River 

HP Hickory Point 

WL Woodlawn 

CB Corbandale 

FCR Fort Campbell Restricted Area 
NC Numerous Creeks 

NSP Numerous Small Ponds 

F Forest and Brushwood 

C Cropland and Pasture 
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The Fort Canpbell restricted zone encompasses over 100,000 acres. 
Except for the restricted area and the forested region, tobacco 
and general tarming dominate the area around Fort Campbell and 
Clarksville. 

A high pressure system centered over northern West Virginia 
produced generally clear skies on 22-23 August 1978. The daytime 
temperatures were around Sl^C and the nighttime temperatures averaged 
about 14”C over the region. No significant precipitation had 
occurred in the area for at least a week. 

The satellite Images were taken twelve hours apart on 22 August 
1978 at approximately 0300 LT and 1400 LT. The model output was 
for 1400 LT on 22 August and 0300 LT on 23 August. This reversal 
of orbit times is justified because it is believed that the 
temperature response between the two days was stationary. 

Table II lists some of the initial parameters used for the 
Land Between the Lakes and Clarksville model runs. The 1200 GMT 
Johnsonvllle, Tennessee soundings for 22 and 23 August 1978 were 
averaged. The average sounding, shown in Figure 18, was used as 
input to the boundary layer model. The average sounding was 
believed to be more representative of the two days. 

4.3 The Land Between the Lakes August 22-23, 1978 

4.3.1 August 22 Daytime Temperatures 

A 12®C variation of daytime surface temperatures is evident in 
Figure 19. The maximum temperature (36°C) is located just west 
of the Blood River (BR) which is only 26®C. This variation in 
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Table II. Initial Model Input Parameters for Land Between the 
Lakes and Clarksville 


Ground Temperature at the 287*K 

Base of the Layer 

Roughness Length 25 cm 

Preclpltable Water 1.88 cm 

Albedo .19 

Emlsslvlty 1.0 

Wind Speed at z ^ 3 m 5.8 msec 

Direction 45“ 



Figure 18. The 1200 GMT Johnsonville, TN average sounding for 
22 and 23 August 1978. 










tamparatur* it in direct rasponta to a change of land uaa, cloaa 
cropped vegetation or urban areas to highly vegetated or well<^atered 
aurfacca. The lower temperatures are recorded over the river end 
Che forested areas (F); whereas, the highest temperatures are found 
over the croplands (C) . Temperatures greater than 29*C are found 
east of the Blood River, south of Jonathan Creek (JC) where there 
is some cropland crossed by a state highway, southeast of the Little 
River (LR) over farmland, and in the central part of the Land Between 
Che Lakes where the town of Golden Pond lies. Golden Pond is a town 

with fewer chan 1,000 Inhabitants, so it's astonishing that a temperature 
maximum is observed. The 32*C isotherm south of the Little River 
is associated with the town of Cadiz and the intersection of two 
major highways. Cadiz has a population of about 2,500 people; 
however, an Intense temperature maximum is not observed here because 
some cloud (cirrus) was located in that region, chough the cloud was 
removed from the analysis by replacing cloud pixel values with 
Chose of neighboring clear sky pixels. The two areas of 30*C north 
of the Little River (LR) conform to patches of cleared areas in the 
forest. The three 32*C isotherms north of the Fort Campbell (FC) 
border are associated with two interstate highways. The three 32*C 
isotherms north of Che Blood River (BR) are closely connected with 
Che areas that have been cleared In Che forest for highways. 

Temperatures ranging from 24*C to 29*C are found over the Land 
Between the Lakes which Is a predominantly forested area. 

The lowest temperatures are found over Che bodies of water 
where minima of 24“ C are found over Kentucky Lake (KL) , Lake 
Barkley (LB), and the numerous small ponds (NSP). 


Si 


4.3.2 August 23 Nighttlmt T«&p«r«tur«s 

The nighttime surface temperatures shown in Figure 20 range 
from 11*C to 23*C. Isotherms valued at 20*C to 23*C are found over 
Kentucky Lake and 18*C to 21*C isotherms are found over Lake Barkley* 
compared with 24*C to 28*C daytime temperatures. Thus, the diurnal 
temperature variation over water surfaces is small. Oke (1978) 
gives the following reasons why water exhibits very little thermal 
response. Water allows short-wave radiation to be transmitted to 
considerable depths and the existence of convection and mass transport 
by the fluid motions permits the absorbed heat to be spread through- 
out a large volume. Also, the unlimited water availability provides 
an efficient latent heat sink, and evaporative cooling tends to 
destabilize the surface layer and further enhance mixing. Finally, 
the thermal capacity of water Is exceptionally large. In fact, 
about three times as much heat Is required to raise a unit voltime 
of water through the same temperature Interval as most soils. 

The nighttime temperatures over the Land Between the Lakes 
range from 13‘C to 14‘C. The croplands exhibit 12*C to 13*C 
temperatures and an 11”C minimum Is observed In the Fort Campbell 
restricted area. 

4.3.3 August 22-23 Moisture Availability 

The Isopleths of moisture availability shown In Figure 21 are 
largest over the dense canopy of trees which predominate the Land 
Between the Lakes region. M approaches 1.0 over a large portion 
of this region. The three minima of .6 located south of Jonathan 



LAND HTWHN 
TNI LAKBS 


5 


Figure 20. The Land Between the Lakes nighttime blackbody surface 

temperature analysis for 0300 LT on 23 August 1978 (®C). 
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Cr««k (JC) «r« •••oci«ttd with a ■t«c« highway and cropland. Tha 
.5 and .6 laoplatha located north of tha Blood Rivar (BR) corraopond 
CO pravloualy foraatad areas that had bean cleared. Values of .7 
are found over tha southwestern portion of the Blood River and M 
IS .4 about two klloaeters west of the river. Another large 
variation of M with land use is found southeast of the Little River 
(LR) where M varies from .8 to .4 over a distance of four kilometers. 
The non-forested, open areas of Fort Campbell also reflect the lower 
moisture availability values. M varies from .6 to .8 in this area. 

Comparing the daytime surface temperatures with moisture 
availability, one can see the :’nverse relationship between theec two 
variables. The highest daytime temperatures correspond to the lowest 
values of H. For example, the 34*C to 36*C isotherms west of Che 
Blood River correspond to tha .4 isoplech of M. Conversely, the 
24 *C to 25*C isotherms over the southern portion of the Land Between 
the Lakes correspond to values of M that approach 1.0. 

4.3.4 August 22-23 Thermal Inertia 

The values of thersutl inertia depicted in Figure 22 range from 
.03 to values which are greater than .07. The elimination of 
values greater than .07 in Che analysis tends to remove detail 
around the lakes where P exceeded .07. The values of P range from 
.05 to .07 over the predominantly forested Land Between the Lakes. 
The .04 minlnuffl south of Jonathan Creek (JC) corresponds to the 
state highway and cropland. The crop and pasture land west of the 
Blood River (BR) and north of Fort Campbell (FD display values of 
.03. Isopleths of .06 are found over the numerous ssiall ponds (NSP) 
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emst of the Little River. The non-foreeted, open arees of Fort 
Campbell display values between .03 and .06. 

The .05 to .06 values of P located over the Land Between the 
Lakes correspond to nighttime temperatures of 13*C and 14"C. The 
highest nighttime temperatures occur over regions vrhere the highest 
thermal inertia values are observed, such as Kentucky Lake and 
Lake Barkley. The largest diurnal temperature change is found 
where the thermal inertia is at a minimum. For example, the .03 
minimum P value over the cropland west of the Blood River corresponds 
to a 23*C diurnal temperature variation. 

4.3.5 August 22 Surface Heat Flux 

The surface heat flux values, shown in Figure 23, are lowest 

over the bodies of water and the forested areas. H varies from 

o 

-2 -2 

35 Wm to 95 Wm over the Land Between the Lakes. The local 
-2 

maximum (95 Wm ) value located over the central portion of the 

Land Between Che Lakes corresponds to the small town of Golden Pond. 

_2 

The local maximum (135 Wm ) south of Jonathan Creek (JC) is 

_2 

associated with a state highway. The local maximum (155 Wm ) 

north of the Blood River (BR) conforms to previously forested 

areas that had been cleared for highways. The largest heat flux 
_2 

values (175 Wm ) are found over the closely cropped vegetation 

west of the Blood River and north of Fort Campbell (FC) . The 
-2 

135 Wm Isopleth south of the Little River (LR) is associated with 

Cadiz. The numerous small ponds (NSP) east of the Little River 

-2 -2 

display values of 35 Wm , with a 15 Wm minimum over one of the 
larger ponds. 
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4.3.6 August 22 Surface Evaporative Heat Fltix 

The isopleths of the surface evaporative flux« shown in Figure 24, 
range from 290 Wm to 610 Urn . The largest values of are found 
over the dense forest canopy where moisture is absorbed through the 
roots of the many trees and then transpired. Budyko (1963) estimates 
the relative value of transpiration to represent 30 to 40 percent 
of the total evaporation even for relatively low vegetative covers, 
such as grass at a height of 10 to 20 cm. The relative value of 
transpiration is presumably even greater from the soil in a forested 
area. Spittlehouse and Black (1979) used Bowen ratio and eddy 
correlation measurements to determine the evapotranspiration of a 
forest. They found evaporation rates corresponding to values of M 
that exceeded .5. 

_2 

Lower evaporative flux values (490 Wm ) are found south of 
Jonathan Creek (JC) over a state highway, south of the Little River 

(LR) over Cadiz, and east of Donaldson Creek (DC) over a highway. 

-2 “2 

values varying from 410 Wm to 450 Mn are found over farmland. 
Evaporative flux minima are found over Kentucky Lake and Lake Barkley 
because of the cool daytime surface temperatures. 

The largest E^ values are found over regions where was at 
a minimum and M was at a maximum. Reginato et al. (1976) experimented 
with Avondale loam and found that when the soil water content is high, 
the hydraulic conducting property of the soil is relatively high, 
thereby allowing enough water to flow to the soil surface to meet 
the meteorological demand of the potential evaporation rate. When 
the soil water diminishes below a certain value (i.e. second stage 
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of drying) » the potential evaporation rate cannot be net, and the 
relative evaporation rate declines with the decreasing water content. 

4. A Clarksville August 22*23, 1978 

4.4.1 August 22 Daytime Temperatures 

The Isotherms depicted in Figure 25 range from 36*C to 23*C. 

The warmest area Is Clarksville where the temperature varies from 
26”C over the southeastern comer of the city to 36‘C over the north- 
western portion. Fort Campbell displays a temperature maximum of 
33*’C. Three very small towns, all with populations of fewer than 
1,000 people, are referred to in the key because, surprisingly, they 
display temperature maxima. Woodlawn (WL) and Corbandale (C6) 
display temperature maxima of 30*C and Hickory Point (HP) displays 
a maximum surface temperature of 33°C. 

Minimum surface temperatures are found over forested areas and 
wetlands (F) . In these regions, the temperature varied from 23°C to 
28°C. The temperature over the numerous small ponds (NSP) and 
numerous creeks (NC) varied from 26®C to 28“C. Cropland and pasture 
areas (C) exhibit temperatures which range from 28®C to 30"C with 
a few ZO’C minima. 

4.4.2 August 23 Nighttime Temperatures 

The nighttime temperatures displayed in Figure 26 show a 5*C 
variation across the area. The maximum temperature (IS^C) is found 
over the numerous small ponds (NSP) and creeks (NC) that are prevalent 
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Figure 25. The Clarksville blackbody daytime surface temperature 
analysis lor 1400 LT on 22 August 1978 (®C). 
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in this part of Tannassaa. Tha only urban araa that shows soma sort 
of haat Island is Clarksvllls vhare tha tamperatura is 14*C. Fort 
Campball does not display alavatad tamparaturas bacausa this is a 
military installation with no factories or industrial areas which 
are important contributors to the nighttime urban heat island. The 
Fort Campbell restricted area (FCR) shows a localized 10*C minimum 
because of the bare ground found here. The forested areas display 
temperatures ranging from ll^C to 14”C, and 11**C to 12*C minima 
are found over closely cropped vegetation where the soil maintains 
a low heat capacity. 

4.4.3 August 22-23 Moisture Availability 

The wide range of moisture availability, found in Figure 27, is 
indicative of varied land use. The highest values of M are found 
over the forested areas (F) and over the numerous small ponds (MSP) 
and creeks (NC) where M approaches 1.0. Clarksville shows a sharp 
variation of M from .9 over the outskirts of the city to .4 over 
the inner portions of the city. Corbandale and Woodlawn show 
minima of .7, and Fort Campbell and Hickory Point both display 
.5 minima. The Fort Campbell restricted area (FCR), a non-forested 
area, displays values ranging from .6 to .9 depending on the amount 
and type of vegetation that were present. High values of M are 
found over the densely vegetated crops northeast of Clarksville; 
however, the values were much lower over the croplands west of Fort 
Campbell. This reflects how M varied with the typa of crop grown 
and the amount of irrigation that was done. 



Figure 27. Ti;e Clarksville moisture availability analysis for 
22-23 August 1978. 
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The negative correlation between daytitoe temperaturea end 
noiature availability is again evident. The .4 minine over Clarksville 
correspond to the 34**C end 36*C aaximua daytine iaothents (Figure 25); 
likewise, the .5 minima over Fort Campbell and Hickory Point (HP) 
correspond to the 33*C maximum isotherms. Conversely, the lowest 
temperatures on the map (23”C) correspond to values of M that approach 
1 . 0 . 


4.4.4 August 22-23 Thermal Inertia 

The thermal inertia values shown in Figure 28 vary from .02 to 
.07. The lowest values of P, surprisingly, are found over the urban 
areas. One would anticipate a maximum of P over the urban areas 
which would reflect the enhanced ability of the urban substrate to 
store and conduct heat. However, minimum values of’ .02 are found 
over Clarksville and Fort Campbell. Low thermal inertia values 
were also found by Dodd (1978) over urban areas. Hickory Point (HP) 
and Woodlawn (WL) display values of .03, and Corbandale (CB) shows 
a higher value of .04. The forests display thermal inertia values 
ranging from .05 to .07. 

Thermal Inertia varied In a similar manner to the nighttime 
temperature pattern. The large values of P (.06-. 07) corresponded 
to the maximum nighttime temperatures (14-lj°C). Lov; values of P 
are found over Clarksville, Fort Campbell, and the other small 
towns; because of the low values of P there, no significant increase 
In nigh .time temperatures is found for these urban areas. 





67 


4.4.5 Auguat 22 Surface Haat Flux 

Figure 29 ahowa a wide variation of cha aurfaca haat flux with 

land uaa. Tha largaat valuaa of are found ovar Clarkavllla 

(170 Wb“^), Fort Caapball (150 Wb’^), Hickory Point (130 Vm~h p 

_2 

and Uoodlawn (110 Wb ). Thaaa araaa corraapond to the aaxlBua 

daytiae iaothanna and ainlauB valuaa of moiatura availability. Haat 

flux Binlota are found ovev tha foraatad araaa where the heat flux 

-2 

values are as low as 10 Un The creeka and ponds also display 
low heat flux values. 

4. 4. 6 August 22 Surface Evaporative Heat Flux 

The values for the evaporative heat flux, shown in Figure 30, 

-2 -2 

range fron 430 Wb in Clarksville to 630 Wb over the forests. 

-2 

A local evaporative heat flux Binimua of 510 Wb is found over 

-2 

Fort Campbell. The 470 Wm isoplech northwest of Fort Campbell 

corresponds with the runways of the airfield located there. Other 

-2 

^ow values of E are f outul over Hickcrv Point (470 Wb ) , Woodlawn 
o 

(550 Wb ^), and Cortandale (550 Wb ^). The largest values of 

-2 

are found over the forests where values of 630 Wm are found. 

-2 

The relatively high values (630 Wm ) found near the Cumberland 
power plant occur because the power plant is situated in an open 
area bordered by trees on three sides and the Cusberland River on 
the open side. Large E^ values are also found over the numeroiis 


creeks, ponds, and croplands. 
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4.5 Cumberland Albedos 

The base map used for the Cundierland albedo study, along with 
key coded areas, is shown in Figure 31. The flight path for 
17 August 1978 of The Pennsylvania State University's airplane is 
superimposed over the map, and this will be referred to in section 4.6. 

Surface albedo values for approximately 1400 LT on 22 August 1978, 
which were derived from visible satellite measurements, are shown in 
Figure 32. The surface albedo is defined as the ratio of the 
reflected solar energy to the Incident solar energy at the earth's 
surface. A simple transfer model for the visible spectrum was used 
to convert the visible satellite DN values to albedos. The satellite 
response in the visible range was 0.56 ym to 1.04 ym. 

The lowest albedos are found over Kentucky Lake and Lake Barkley 
where values from 4% to 8% are found. Smoothing of the data removed 
lower albedo values over these lakes. The forested areas exhibit 
values from 10% to 14%. The spatial fluctuations of albedos over 
the forest are due to the variability of the forest density, the 
woodlot spacing, and the mixing of tree species. The farmlands 
north of Fort Campbell reveal 14% to 16% albedo values depending 
on the type of crop being grown. The non-forested Fort Campbell 
restricted zone (FCR) displays albedo values ranging from 10% to 
16%. 

These satellite derived albedos agree well with typical 
albedos for various surfaces which are given by Sellers (1965) 
and Kung et al. (1964). Rung et al. (1964) found the albedo values 
of forests to be lower than that of farmlands. The following typical 
albedos are given by Sellers (1965) for various terrain: deciduous 




Figure 31. The base map used for the Cumberland albedo study and 
the flight path for 17 August 1978 of The Pennsylvania 
State University’s airplane, along with the following 
key areas coded for easy reference: 

LB Lake Barkley 

KL Kentucky Lake 

CR Cuinberland River 

CB Corbandale 

WL Woodlawn 

FCR Fort Campbell Restricted Area 
NSP Numerous Small Ponds 

NC Numerous Creeks 

F Forest and Brushwood 

C Cropland and Pasture 
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' Figure 32. The Cumberland surface albedo analysis (derived 

from visible satellite measurements) for 1400 LT 
on 22 August 1978 {%) . 
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forests, 10— 20X; coniferous forests, 5— 15X; crops, 15 — 2SX; green 
meadows, 10— ‘20X; black topped roads, 5— lOX; dark soil, 5— 15%; 
and water, 3— 5X at solar noon. 

Rung et al. (1964) observed that the surface albedo is related 
to the phenological cycle of the vegetation. Lush tree leaves and 
lower vegetation tend to refle t a greater fraction of incident 
solar energy than do canopies with sparse tree leaves or vegetation. 
Soil moisture content may also be a probable factor Influencing the 
surface albedo. Incident rays can enter in all directions when the 
soil particles and plant surface are covered with a film of water; 
however, only rays which can reach the surface of the water film 
within the limiting angle of total reflection can emerge. The 
surface albedo can also be directly influenced by the soil color 
and water holding capacity; thus, soil types should be considered 
when looking at regional surface albedos. However, Rung et al. (1964) 
believe that the surface albedo difference due to soil types is less 
important than that due to vegetative covers, at least where bare 
soils are not comnon. The surface albedo is also dependent on wave- 
length. The wavelengths ranging from 0.5 ym to 1.0 ym (near infra- 
red) depend on the solar angle and the viewing angle of the satellite. 

4.6 Comparison of Aircraft Measurements with Satellite Derived 

Albedos, Surface Heat and Evaporative Fluxes 

From 10 August to 19 August 1978, The Pennsylvania State 
University's Aerocommander 680E made nine flights for the STATE 
project. A detailed description of these flights is given by 
DiCristofaro and Albrecht (1979). The flight path for one of these 
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flights is displayed in Figure 31. This flight was a land use study 
made in the vicinity of the Cumberland power plant. The passes were 
made at 500 feet (MSL). 

Eight radiometers on board the aircraft were used to measure 
the short and long-wave radiation. The ratio of the reflected short- 
wave Irradlances to the incident irradlances was used to calculate 
the albedo of the surface. Simpson's Rule (described in the Appendix) 
was used to calculate the average albedo for each flight pass. The 
surface albedos derived from the visible satellite measurements, 
shown in Figure 32, were also averaged for each flight pass through 
the use of Simpson's Rule. 

A comparison of the aircraft and satellite measured albedos for 
each flight path is given in Table III. The aircraft measured values 
were 1.5% higher on the average; thus, the agreement was excellent. 

The largest differences were for pass number three (2.7% difference) 
and pass number six (2.4% difference). These were the shortest of 
the nine passes. The best agreement (0.4% difference) was for 
pass number five which was one of the longest passes. Better 
results are obtained if a bidirectional transfer model for the 
visible spectrum is used to convert the visible satellite DN values 
to albedos. 

From 2 August to 28 August 1978, the Boundary Layer Group (BLDG) 
of Che National Oceanic and Atmospheric Administration (NOAA) made 
airborne turbulence measurements in the vicinity of the Cumberland 
power plant. Four of the flights (9, 10, 14, and 27 August 1978) 
were land use studies made at an altitude of approximately 250 to 


300 m. 
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Table III. The Satellite Measured Albedos for 22 August 1978 vs. 

The Aircraft Measured Albedos for 17 August 1978 at 
Approximately 1400 LT. 


Pass 

Ntunber 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Satellite 
Albedo (X) 

12.4 

11.5 
11.0 
11.8 

12.5 
11.7 

11.3 

11.4 
11.7 
12.1 


Aircraft 
Albedo (X) 

13.4 

13.2 

13.7 

12.7 
12.9 
13.1 

13.7 

13.1 
13.4 

13.2 


A comparison of the aircraft measurements for the four different 

days with the sacelllce measurements for 22 August 1978 was made. 

The aircraft measurements were corrected to surface values by 

multiplying the aircraft values by where h is the height of 

a 

the mixed layer and is the height of the aircraft. A correction 
for the time difference was also made. Ching (1980) maintains that 
it is not valid to neglect synoptic variations; therefore, a 
comparison of the satellite measurements for 22 August 1978 should 
not be made with the aircraft measurements for the four different 
days. Indeed, the largest differences between the aircraft and 
satellite measurements were for 9, 10, and 14 August 1978. The 
closest agreement was for 27 August 1978 and these results are 
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given in Table IV. The filtered aircraft results were obtained from 

NOAA. The satellite values were obtained from Figures 23 and 24, 

and Simpson's Rule was used to obtain an average value of H and E 

o o 

over each flight path. The passes were made over the Land Between 
the Lakes. 

Table IV. The Satellite Measured Heat and Evaporative Fluxes for 
22 August 1978 vs. The Aircraft Measurements for 
27 August 1978. 


Pass 


Satellite 

Aircraft 

Number 


Value (Wm~^) 

Value (Wm‘ 

■) 

Eo 

498.5 

450.5 

L 

Ho 

68.0 

145.4 

1 

Eo 

527.2 

489.7 

0 

o 

67.0 

198.7 

5 

E 

533.2 

512.0 

o 

65.2 

163.8 


The aircraft measured values of the evaporative heat flux 

-2 

were, on the average, 36 Win lower than the satellite measured 
values. The large heat flux values obtained by the aircraft 
measurements seera too large and the spatial varir* i t i^at flux 
is opposite to what is expected. The satellite images shown in 
Figures 3 through 6 and the daytime surface temperature analysis 
shown in Figure 19 clearly display lower temperatures over the 
Land Between the Lakes area; yet, the aircraft measured heat flux 
values were found to be larger in this region. The aircraft heat 
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fluxes over the forested areas were approximately the values found 
for downtown St. Louis (see Figure 14)! Dense vegetation, such as 
the canopy of trees found In the Land Between the Lakes region, is 
known to contribute significantly more evaporative than sensible 
heat flux, l.e. the Bowen ratios are typically small. 

4.7 Summary of Results 

A large variation of the daytime surface temperature was 
observed over the Land Between the Lakes and Clarksville areas. 

Small towns with populations less than 1,000 people exhibited local 
temperature maxima. The lowest daytime temperatures were observed 
over the forested areas where a large amount of the evaporation was 
also observed, thus causing the cooler conditions. The highest 
daytime temperatures corresponded to the areas of lowest moisture 
availability and evaporative fluxes, and highest surface heat 
fluxes. Evaporation is suppressed on dry surfaces whlca allows a 
greater partitioning of the net radiation into sensible heat, thus 
causing greater surface heating. More energy is partitioned Into 
latent heat flux for moist surfaces causing suppressed heating. 

The largest diurnal temperature range was found over farmland 
where the thermal inertia was a minimum. The smallest diurnal 
temperature variation was found ever the bodies of water. Maximum 
values of thermal inertia were anticipated over the cities of 
Clarksville and Fort Campbell. Tnis would have reflected the 
enhanced ability of the urban substrate to store and conduct heat. 
However, minimum values of thermal inertia were found over some 
urban areas such as Clarksville and Fort Campbell; thus, no distinct 
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nighttime temperature maxima were evident. It is evident that 
warmer nighttime temperatures are also related to higher temperatures 
during the day, not just to P. 

The distinct variation of vegetation over the area, or a 

variation of transpiring surfaces, was responsible for a large 

change in the heat and evaporative fluxes. Large changes of and 

were observed over a distance of even a few kilometers. Tnus, 

vegetation is the arbiter of H and E . 

o o 

The moisture availability values a.-proached 1.0 over the 
forested areas which accounts for the large surface evaporative 
fluxes found there. Shuttleworth and Calder (1979) made evaporation 
measurements for tall (forest) vegetation for two cases. The 
comparison was made for conditions when water was freely available 
on the surface of the vegetation, and for a "potential" condition 
when there was no significant moisture stress, which is probably 
an erroneous condition. Moisture availability values of .5 were 
found for the dry conditions and values of 1.0 were found for the 
wet conditions. Shuttleworth and Calder (1979) concluded that the 
status of the surface (i.e., whether it is wet or dry) within a 
given type of vegetation, is the primary control on the evaporation 
rate. 

The albedo values over the Cumberland area were lowest over 
the bodies of water and highest over the croplands. The albedos 
varied by 6% over the various vegetative covers. Agreement between 
aircraft derived values and those determined from the satellite 
measurements was close, although there was very little spatial 


variation in the former. 
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5.0 INFERRING THE EFFECTS OF SURFACE HEAT FLUX ON 
PLUME SPREAD AND CONCENTRATION OF CONTAMINANTS 


5.1 Atmospheric Dispersion Resesrch 

Presently, the most popular method of determining a and a is 

y * 

from the Pasquill-Gifford empirical curves in which 0^ and 0^ are 
plotted as a function of x. As Panofsky (1969) points out, the main 
drawback of the Pasquill-Gifford curves is that they do not allow for 
the effect of terrain roughness. The empirical curves were based 
on experiments over smooth terrain; therefore, the dispersion over 
cities and other rough regions is underestimated. The Pasquill- 
Gifford curves were produced from ground sources and represent 
averages over a few minutes. Thus, the curves must be adjusted if 
they are to be used under conditions different from those for 
which they were constructed. 

The classic model used by most atmospheric dispersion scientists 
for elevated point sources is the Gaussian plume diffusion model. 

This is described by 


2 2 
X - ^ exp (- -^) (exp (- 

2TTOyO^U 20y 


(5) 


, (z+H)'", , 

+ exp ( =-) ] 

20 ^ 
z 


where x is the ground-level concentration, Q is the source strength, 
0 ^ and a are the standard deviations of the time-averaged plume 

V* 2 

concentration distribution in the y (cross-wind) and z (vertical) 
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directions, U is the mean horizontal wind, and H is tha aff active 
source enission height. 

Tha Gaussian plume model is based on steady-state conditions, 
homogeneous flow, no atmospheric chemical reactions or gravity 
fallout, a sampling time of at least ten minutes, and a significant 
wind speed. Also, for the form of the Gaussian model presented in 
equation (5) , perfect reflection of the plume at the ground is 
assumed . 

Munn et al. (1972) point out that equation (S) does not apply 
when the winds are light, when serious air quality deterioration 
can occur. The Gaussian plume dispersion swdel also assumes that 
there are no vertical gradients In wind. Also, steady-state 
homogeneous conditions may be a mathematical convenience rather 
than reality. The Gaussian plume model does not directly take into 
account a variable heat flux. In some dispersion studies, such as 
the one conducted by Venkatram (1980), the surface heat flux was 
assumed to be 0.3 times the incoming solar radiation. However, it 
will be shown in section 5.3 that a variable heat flux has an 
important effi^ct on plume spread and concentration. 

Many new methods of estimating dispersion parameters have been 
suggested in the past few years. One such method was developed 
by Willis and Deardorff (1974). They studied the turbulence and 
transfer properties in the laboratory under free convection 
conditions. Free convection occurs when turbulence is independent 
of surface drag force such as in the daytime under cleat skies and 
light winds when dispersion is often dominated by convection. The 
surface heat flux enters the Willis and Deardorff method explicitly; 
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therefore, chie method ie ueed rether then the Peequill-Glfford 
method. The Willis and Deardorff method of estimating dispersion 
parameters is discussed in detail in section 5.2. 

5.2 Dispersion from an Elevated Point Source in the Convective 

Planetary Boundary Layer 

Willis and Deardorff (1978) used a laboratory model and 
Lamb (1978) used a numerical model of the convective planetary 
boundary layer to study dispersion from an elevated source. The 
results applied to a continuous point source located within a 
thermally convecting field in the presence of a simulated uniform 
mean wind. Free-convection similarity scaling was used to non- 
dimenaionalize the quantities measured. An elevated release height 
of * .25 h was chosen by Willis and Deardorff (1978) and 
z^ ■ .26 h was chosen by Lamb (1978), where h is the thickness of 
the convectlvely mixed layer. 

The simple, workable laboratory model consisted of a deep 
convective layer heated from below with a capping inversion formed 
by an overlying stable layer. A detailed description of the model 
is given by Willis and Deardorff (197A). 

The turbulent velocities and in the upper part of the 
planetary boundary layer scale with the convective velocity w^ given 
by 


w 




'oc t' 
P 


9 


( 6 ) 
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where is the surface heat flux, p is Che air density, is the 
specific heat at constant pressure, and T is the average temperature 
of the mixed layer. The controlling scales in the mixed layer are 
w^ and h. This scaling is valid provided the height exceeds L, 
given by 



(7) 


where k is Che von Karman constant and u^ is Che friction velocity. 

These scales are then used to nondimensionalize all the 
variables. The dimensionless time xised is 


* 


c 



( 8 ) 


where t is the time measured from the centered time of particle 
release. Taylor's hypothesis, t *>’ , is used to convert time to 

a space coordinate, where x is the distance from the pollutant 
source and U is Che simulated uniform mean wind blowing along Che 
X axis. The dimensionless time or distance can Chen be written as 

t -x-(— )(^) . (9) 


This transformation Is valid within the limits of 


1.2 w^ < U < 6 w^ 


( 10 ) 
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Willis and Desrdor££ (1978) maintain that tha lowar limit ia naedad 
to insura that tha maan wind is strong anough so that signi£icant 
upstraam turbulant di££usion doas not occur. Tha uppar limit 
represents the value above which turbulence in tha bulk o£ tha mixed 
layer begins to be signi£icantly a££ected by shear stress. 

The nondimensional mean particle height is de£ined as 

fl.lh 

Z - 2 p dz . (11) 

^0 * 


where p^ dz is the probability that a particle resides in a 
horizontal layer between z - -|Az and z -|az. The numerical results 
of Lamb (1978) are shown in Figure 33 where the nondimens ionalized 
mean height Z is plotted as a £unction of the nondlmensionalized 
downwind distance from the source X. Initially, Z remains 
approximately constant at short distances of X until enough 
descending particles have been entrained into updraft regions. The 
Gaussian plume or K-Theory models of turbulent diffusion cannot 
describe this aspect of diffusion in a convectlvely mixed layer. 

Z then increases to its equilibrium value of about 0.5 h at X « 1.6, 
with a slight overshoot of approximately 0.025 h. 

The dimensionless vertical and lateral, overall standard 
deviations of particle displacement are defined, respectively, as 


(Z 


- a . [ 


rl.lh 


,, , .2 . , 1/2 

(Z-Zg) p^ dz] 


( 12 ) 


(Y 


. a - ( 


r2h 


v2 . ,1/2 

Y Py dy 1 


-2h 


(13) 
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Figure 33. The nondinenslonal mean particle height (Z) as a 
function of the nondlmenslonal downwind distance 
from the source (X) (from Lamb, 1978). 



Figure 34. The dimensionless overall standard deviations as 
a function of the dimensionless downwind distance 
from the source (from Willis and Deardorff, 1978). 
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wh«r« Py dy is tht probability that a partlcla vlll appear anywhara 
in a vartical column batvaan y - <|i6y and y ^ -l^y. Tha laboratory 
raaulta» mada by Wlllia and Daardorff (1978), of tha dimanaionlaaa 
overall particle atandard daviationa are ahown in Figure 34 aa a 
function of the dimanaionleaa downstream distance. At tha source 
height, the root-Bean~8quara turbulent vartical velocity, 0^, is 
greater than the turbulent lateral velocity component, therefore, 
the vertical spread, a , exceeds the lateral spread, o , at short 

* y 

values of X. The effects of the capping inversion begin to slow 
the vertical spread by X > 0.8, until reaches its equilibrium 
value of about .38 h with only a slight overshoot. The lateral 
spread of the plume continually Increases with X, with a sharp 
increase near X - 2.2. Willis and Deardorff (1978) believe this 
feature is due to some of the particles being trapf«:d for a time in 
Che lateral confines of a large-scale eddy befove undergoing 
accelerated lateral spread once sufficient numbers of particles are 
carried into an adjacent large-scale eddy. 

The dimensionless mean cross-wind integrated concentration is 
defined as 


F (X,Z) 


P' (X.Z) hu 

Q 


(14) 


where 




(X.Z) 


lim 

t-*«o 


p 


c (X,Y,Z,t) dY 


(15) 
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The time dtptndcnt point Man concentration is c(X,Y,Z,t) and the 

aaiasion rate is Q. The x, y, and z coordinates nondiMnsionalized 

by h are represented by X, Y, and Z in equations (14) and (IS). 

The nusMrical simulation of ^ perfomted by Land) (1978) is shotm 

in Figure 35. This is for an elevated source height of z * 0.26 h 

and a time average over a period of about 20 minutes or longer. The 

centerline of the plume descends from the source to a point near 

X 'V, , where the centerline intercepts the ground, until a 

distance x ^ . The plume then begins to ascend again. 

Willis and Deardorff (1976) computed the cross-wind integrated 

concentrations for a source height of z ■ 0.0h7 h. After comparing 

8 

the results of the lower release height with the results computed 

for z *0.25 h, they concluded that the source height is a dominant 
s 

influence until the concentration maximum first reaches the ground. 
Beyond that distance, the concentration fields are very similar. 

The ground- level values of the mean cross-wind integrated 
concentration, ^(X,0) , as computed by the numerical study made by 
Lamb (1978), are presented in Figure 36. Comparing these results 
with the Gaussian plume results. Lamb (1978) found the ground- 
level concentrations to be larger than chose predicted by the 
Gaussian plume formula at all distances from the source. 

5.3 St. Louis and Cumberland Plume Characteristics 

The heat flux maps for St. Louis, the Land Between the Lakes, 
and Clarksville shown in Figures 14, 23, and 29 display a large 
variation of heat flux over a distance of a few kilometers. The 
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a 2 a4 as ae l 

I 
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Figure 35. The dinenelonless ncan cross-wind integrated 

concentration (^) for an elevated source height 
of Zg"0.26h indicated by the vertical bar. The 
centerline of the plume is represented by the 
dashed line (from Lamb, 1978). 



Figure 36. The ground-level values of the mean cross-wind 
integrated concentration (^(x,0)) from a 
continuous point source of height 2 ^* 0 . 26h as 
a function of the dimensionless downwind 
distance from the source (from Lamb, 1978). 



j 

I 

I 



89 


purpose of this study is to invsstigsts ths inpset of this vsrisbls 

host flux on ths pluas sprssd and concentration. 

Typical values for the height of the aixed layer (h), wind 

speed (U), aean tenperature of the layer (T), density (p)» and 

specific heat (c ) used to compute and X are listed in Table V. 

P " 

The effect of a variable heat flux on plume spread and concentration 
is studied at a distance of 1 km and 2 km from the source. The 
results for x ■ 1 km are given in Table VI and the results for 
X ■ 2 km are given in Table VII. The values of heat flux listed 
in Tables VI and VII are representative values found over St. Louis 
and Cumberland during the period studied. 

Table V. The Values Used to Compute and X. 

h ■ 1000 m 
U - 2 mser"^ 
p ■ 1.2 kgm ^ 

T • 299’K 

Cp • 1004 Jkg”^ K"^ 

The following discussion is in reference to the results 
presented in Tables VI and VII. These results were obtained by 
the method described in section 5.2. 

The nondimens ional mean particle height (Z), at x " 1 km, 

-2 -2 
increases from .28 for H ■ 10 Wm to .38 for H ■ 200 Wm . 

0 o 

The mean particle height, at x ■ 2 km. Increases from .34 for 

H - 10 Wm"^ to .51 at H • 200 Wm’^. 

0 o 


Table VI. The Effect of Surface Heat Flux on Various Pluae Spread and 
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2 1/2 

The dlMMionlass vertical ((Z' ) ' or 0^) and lateral 
2 1/2 

((Y* ) or Oy)t overall atandard daviatlona of particle diaplaceaent 
alao increase with Increased heat flux. Under low heat flux 
conditions and at short values of X, o exceeds a at both x ■ 1 l« 

* y 

»2 

and X ■ 2 ka. At x ■ 1 ka« 0^ equals at ■ 100 Wb . The 

capping inversion then begins to slow the vertical spread of the 

-2 

pltine so that o exceeds G at H ■ 200 Wa and x ■ 1 ka. At 
y z o 

-2 

x * 2 ka, G exceeds g at H SO Wa . 

' y z o 

The nondiaensionalized aean cross-wind integrated concentration 
(^) peaks at 1.9 for “ 10 Wa~^. ^ decreases to 0.9 under the 

-2 -V 

largest heat flux conditions (200 Wm ) . At x - 2 ka. varies 

froB 1.0 to .95 with no discernable pattern. 

The ground-level value of the mean cross-wind integrated 

concentration (^(x.O)), at x - 1 ka. peaks at 2.1. This corresponds 
-2 —V 

with a 100 Wm heat flux. At x • 2 km, C‘^(x,0) peaks at a lower 
-2 

hint flux (10 Wm ). The lowest ground- level values, at x ■ 2 ka, 
are found with the highest heat fluxes. 

Lamb (1978) found that the maximum ground-level cross-wind 
integrated concentrations occurred at 


0.6hU 


w. 


( 16 ) 


Using this equation, the maximua value of x (1846 a) was found 

8 

-2 

under the lowest heat flux conditions (10 Wm ). Similarly, the 


minimum value of x (678 a) was found under the highest heat flux 
8 

-2 

conditions (200 Wm ) . 
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Th« Man croaa-wlnd intagratad concantratlona vara \iaad to 
cooq>uta tha Man ground-laval point concantratlona by tha follovlng 
Mthod. Tha Man croaa-wind Intagratad concantratlon la dafinad aa 


CIC 


c 


X dy 


- § F(*) 


(17) 


vhara x 1> tba point concantri^tion, Q ia tha amiaaion rata, and F(x) 
is a dispersion function which dascribas tha growth of tha pluM in 
the z direction. The point concentration is defined as 


X • ^ G(y) F(*) 


(18) 


tdtere G(y) is a dispersion function which describes the growth of 
the pluM in the y direction. Since the cross-wind integrated 
concentration is already known, then one only needs to know G(y) 
to compute the point concentration. G(y) is defined as 


G(y) 




exp(- 


:) 


(19) 


2o 


which reduces tu 


G(y) 




( 20 ) 


for centerline pluM concentrations or at y ■ 0. The Man grotmd 
level point concentration is then 
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X(x,0,0) • . (21) 

L«^ (1978) txpr«ss«d eh« Man point eoneontration in taxaa 
of Q/hU; howavar* ehia la a conatant factor In thla atudy. It la 
ballavad that although tha haat flux changaa draaMtlcally ovar a 
dlatanca of a faw kllooMtara, tha halght of tha ailxad layar la 
dapandant on tha ovarall Intagrated haat flux, ao h can ranaln 

conatant (Panofaky, 1980). 

— -2 
At X • 1 ko, x(x»0,0) decraaaoa from 3.91 at ■ 10 Ua 

•2 — 

to 2.61 at ■ 200 Wa . At x ■ 2 km, x(x«0,0) dacraaaaa from 

3.49 at H - 10 Wm'^ to .91 at H - 200 Wm"^. 
o o 

Figure 37 dlaplays tha maan ground-laval point concentration 
ae a fmctlon of tha dlatanca from tha aource for three conatant 
heat flux valuea. Aa the heat flux Inereaaaa, the cloaer the 
point of maximum ground concentration la to the aource. The 
maximum ground- level point concentration la apparently independent 
of the heat flux and occura at 4.1; thua, expreaaed aa 

W ■ • <”) 

The Inverae heat flux plotted aa a function of the diatance 
from the aource of the maximum point concentration (x ) la ahovn 

BAX 

in Figure 38. Aa (H )”^ 0 and x •►0. Llkewiae, the 

O 0 BAX 

amaller the heat flux, the farther the maxlmiai ground-level point 


concentration ia from the aource. 



0 .S 1.0 13 2j0 23 33 3.9 4.0 43 9.0 9.9 CO 

I Ckfii3 


Figure 37. The mean ground-level point concentration 
(X(xt0,0)) aa a ftmetien of the downatream 
distance from the aource for H.«200, 50, ...>d 

10 Wm’2. 
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S.4 Suasury of Results 

The Willis sad Dssrdorff method of free convection seeling wes 
used to infer the effects of surfece heet flux on plums spread 
end concentration. 

An increased heat flux (H ) causes an increasad convective 

o 

vertical velocity (w^) which was found to cause an increase of the 
nondimensional mean particle height (Z) . The largest increases 
occurred 2 km from the source where the particles have had a longer 
time to rise. 

The increased vertical velocity caused by the increased surface 

heat flux, also caused an increase of the vertical spread (a ) of 

z 

the plume. An increased heat flux was also found to cause an 

Increase of the lateral spread (o^) of the plume. The largest 

values of a and o were found at x ■ 2 km under the largest heat 
y z 

_2 

flux conditions (200 Wm ) . 

The largest mean cross-wind integrated concentrations (^) 

-2 

were found under the lowest heat flux conditions (10 Wm ) where 

there was the least amount of vertical mixing of the plume. The 

largest variation of was found at x ■ 1 km. Farther from the 

source, at x ■ 2 km, the heat flux was not an important consideration. 

At x * 1 km, the largest ground-level values of the mean 

cross-wind integrated concentration (^(x,0)) were found at 
-2 

• 100 Wm . Farther from the source, at x ■ 2 km, the largest 
concentrations were found under the lowest heat flux conditions 
(10 Wm"^). 

The smaller vertical velocities associated with lower heat 
fluxes caused the location of the maximum ground concentration (x ) 


98 


to be farthest from the source. The pollutant was brought down 

closer to the source under higher heat flux conditions. A difference 

in X of over 1 km was found for heat fluxes that varied from 
8 

10 Wm~^ to 200 Wm"^. 

The ground- level point concentrations (x(x,0,0)) were largest 

for smallest heat fluxes due to less vertical mixing. At x ■ 2 kn, 

— -2 
a very large difference of x vas found between 10 Wm and 

200 Wm"^. 

The maximum ground-level point concentrations were 

found to be Independent of heat flux, but not independent of the 
distance from the source. As the heat flux increased, the vertical 
spread of the plume increased so that the point of maximum concentra- 
tion was found closer to the source. For an infinitely large heat 
flux, would be found at the source. 

It has been shown in preceding chapters that heat flux varies 
substantially over different terrain. This chapter has proved that 
this variable heat flux has an Important effect on plume spread and 
concentration. This variable heat flux should be an Important 
consideration when one models dispersion from a power plant. The 

variable heat flux should also be considered when one plans to 

-2 

build a power plant. Large surface heat fluxes (200 Wm ) are 

found over urban and industrial centers; whereas, low surface heat 
-2 

fluxes (50 Wm ) are found over forested areas. The higher the 
heat flux, the closer the maximum ground concentrations are to the 
source. It is therefore important to consider the land use in the 
downwind direction of a power plant. 
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6.0 CONCLUSIONS 


A method has been presented that couples a boundary layer 
model with satellite derived surface temperatures to deduce values 
of omisture availability, thermal Inertia, heat and evaporative 
fluxes. A sophisticated satellite image display system was used 
to determine these various parameters from remotely sensed data. 

This system saved a treoiendous amount of time compared to earlier 
methods of data processing, and Increased the accuracy of the 
final product. 

For the three cases studied, the maximum daytime surface 
temperatures corresponded to the minimum values of moisture 
availability and this was found in the urban centers. The maximum 
nighttime surface temperatures were found to be related to the 
maximum thermal inertia values and this was found over the bodies 
of water and the urban center of St. Louis. Clarksville and 
Fort Campbell displayed minimum values of P; therefore, no distinct 
nighttime temperature maxima were observed. Thus, it was seen 
chat M and P were very much responsible for the formation of 
important temperature variations over Che urban- rural complex. 

The warmer nighttime temperatures were also related to higher 
temperatures during the day. 

Direct measurements of heat and evaporative fluxes are very 
costly and quite cumbersome. A reliable iMthod of measuring 
and remotely without the need for surface or aircraft instrumenta- 
tion has been presented. Maxima of H and minima of E were located 

^ 0 0 


In the urban centers end these regions corresponded closely to the 
sites of msxlauiB deytisM temperatures. The opposite held true for 
the rural erees. Smell Bowen ratios were found over vegetated 
areas and Bowen ratios of approximately 1.0 were found over the 
urban centers. Hence, a large fraction of the land uae variability 
in heat and evaporative fluxes Is caused by variability In rates of 
evapotransplration. 

Satellite determined measurements of heat and evaporative fluxes 
and visible channel albedos were compared to corresponding values 
made by aircraft. In the case of surface heat flux the agreement 
was poor; with surface evaporative flux there was some agreement, 
while the albedos compared very closely. 

Heat flux values are an Important Input to mesoscale and 
diffusion models; however, their values are often surmised. The 
Willis and Deardorff method of estimating dispersion parameters 
was used to infer the effects of surface heat flux on plume spread 
and concentration. Most importantly, it was shown that the ground- 
level point concentrations were largest for the smallest heat 
fluxes; however, the maximum ground- level point concentrations were 
found to be independent of heat flux. independent 

of the distance from the source, though. This variable heat flux 
effect on plume spread and concentration should be an Important 
input to diffusion models. This should also be an important 
consideration in choosing the sites of future power plants. 

Further research is needed to verify the satellite determined 
surface heat fluxes with simultaneous direct measurements. Also, a 


coiq>«rlson of th« disportion paraaotort darlvad by tha Uillla and 
Daardorff aathod with raal diffualon axparlnanta la naadad whan 
tha aurfaca haat flux is known. 


APPENDIX 


DESCRIPTION OF SIMPSON'S RULE 


The satellite derived and aircraft measured albedos listed in 
Table III and the satellite derived heat and evaporative flux 
values listed in Table IV were calculated using Simpson's Rule. 
This can be expressed as: 


" f(x)dx - -| (fg 4fj^ + 2f^ + Afj + 2f^ + 


+ 2f o + 4f , + f ) 
n-2 n-1 n 


where: I • , x > x 

n no 

n * number of subintervals ■ 2, 4, 6, 8, ... 
This is presented graphically In Figure 39. 
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